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ABSTRACT 

We assess the possibility to detect and characterize the physical state of the missing baryons at low redshift 
by analyzing the X-ray absorption spectra of the Gamma Ray Burst [GRB] afterglows, measured by a micro 
calorimeters-based detector with 3 eV resolution and 1000 cm^ effective area and capable of fast repointing, 
similar to that on board of the recently proposed X-ray satellites EDGE and XENIA. For this purpose we have 
analyzed mock absorption spectra extracted from different hydrodynamical simulations used to model the 
properties of the Warm Hot Intergalactic Medium [WHIM]. These models predict the correct abundance of 
OVI absorption lines observed in UV and satisfy current X-ray constraints. According to these models space 
missions like EDGE and XENIA should be able to detect ~ 60 WHIM absorbers per year through the OVII 
line. About 45 % of these have at least two more detectable lines in addition to OVII that can be used to 
determine the density and the temperature of the gas. Systematic errors in the estimates of the gas density 
and temperature can be corrected for in a robust, largely model-independent fashion. The analysis of the GRB 
absorption spectra collected in three years would also allow to measure the cosmic mass density of the WHIM 
with ~ 15 % accuracy, although this estimate depends on the WHIM model. Our results suggest that GRBs 
represent a valid, if not preferable, alternative to Active Galactic Nuclei to study the WHIM in absorption. The 
analysis of the absorption spectra nicely complements the study of the WHIM in emission that the spectrometer 
proposed for EDGE and XENIA would be able to carry out thanks to its high sensitivity and large field of view. 

Subject headings: 



1. INTRODUCTION 

According to the standard cosmological scenario, baryonic 
matter provides a small fraction, about 4-5 %, of the total 
cosmic energy content. While baryons are to first order dy- 
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namically negligible in shaping the large-scale structure of the 
Universe, they have been the essential component in creating 
stars, galaxies, planets and the life in the Universe and consti- 
tute the only component that interacts with electro-magnetic 
radiation. As such, they allow to study in great detail, through 
diverse observations, the physical conditions for the forma- 
tion and evolution of cosmic structures. 

As a matter of fact the baryon density in the Ly-a forest at 
redshift 2-3 (Weinberg et al. 1997, Rauch 1998) accounts 
for a large fraction of the baryon mass as inferred from the 
predictions of the Big Bang Nucleosynthesis combined with 
the measured abundances of the light elements (Buries and 
Tytler 1997, 1998) and, independently, with the temperature 
fluctuations of the cosmic microwave background (Spergel et 
al. 2003, Komatsu et al. 2008). On the other hand, an ob- 
servational census of cosmic baryons in the local Universe 
(Fukugita, Hogan & Peebles 1998, Fukugita & Peebles 2004, 
Danforth & Shull 2005) shows that at z only about 50 
% of the baryons have already been detected. The content 
of stars in galaxies and the interstellar medium provides only 
about 8 % of the overall baryonic budget; the hot plasma with 
temperature of at least 10^ K contained in the potential wells 
of galaxy clusters adds an additional 11%. Finally, the diffuse 
gas with temperature below 10^ K, that constitutes the local 
Ly-a forest, provides a further 29 % (Bregman 2007). 

Hydrodynamical numerical simulations have provided a so- 
lution to this paradox within the framework of the 'concor- 
dance' ACDM model. They predict that about half of the to- 
tal baryonic content of the Universe at z < 1 should be in a 
warm-hot phase with temperature in the range 10^ - 10^ K, 
distributed in a web of tenuous filaments (Cen & Ostriker 
1999a). As such, the WHIM represents the dominant com- 
ponent of baryons in the low redshift Universe. Its thermal 
evolution is primarily driven by shock heating from gravita- 
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tional perturbations breaking onto unvirialized cosmic struc- 
tures, such as large scale filaments and virialized structures, 
such as galaxies, groups and clusters of galaxies. Feedback 
effects related to the nuclear activity or to galactic superwinds 
from galaxy and star formation also play a role, especially in 
high density environments (e.g. Cen & Ostriker 2006). These 
baryons are so hot and so highly ionized that they can only 
emit or absorb in the far UV and soft However, being the typ- 
ical density of the WHIM lower than that of the hot Intra- 
Cluster Medium it is extremely hard to detect. Yet, it must 
contribute to the soft X-Ray Background. The importance of 
this contribution clearly depends on the physical conditions 
of the WHIM (e.g. Voit & Bryan 2001, Bryan & Voit 2001). 
The exquisite spatial resolution of the Chandra X-ray obser- 
vatory has now allowed to resolve a large part, more than 
80 %, of the cosmic X-ray background at soft energies into 
the contribution of point-like sources (mostly Active Galac- 
tic Nuclei - AGN hereafter). Hickox & Markevitch (2007a, 
2007b) recently analysed the Chandra Deep Fields North and 
South. After removing the contribution from X-ray, optical 
and IR sources and carefully characterizing the instrumen- 
tal background and non-cosmological foreground, they com- 
puted the residual background. According to their analysis, 
this residual background in the 0.65-1 keV band cannot be 
accounted for by extrapolating the source number counts to 
lower fluxes. Instead, it comes interestingly close to the theo- 
retical predictions of WHIM emission based on cosmological 
hydrodynamical simulations that include realistic description 
of star formation and SN feedback (e.g. Roncarelli et al. 
2006). If due to soft X-ray emission from the diffuse gas, this 
residual X-ray background would represent the first, though 
indirect, detection of the missing baryons in emission. 

Clearly, the study of the unresolved soft X-ray background 
just provides an upper limit to the amount of diffuse gas and 
indirect indications of its physical properties. The best chance 
to detect the WHIM is by observing the emission or absorp- 
tion lines of highly ionized element like C, N, O, Ne, and 
possibly Mg and Fe. The strongest lines expected correspond 
to the outer K (X-ray band) and L (Far UV band) shell transi- 
tions from H-like, He-like and Li-like Oxygen ions: the OVIII 
ls-2p X-ray doublet [from now on simply OVIII], the OVII 
ls-2p X-ray resonance line [OVII, hereafter] and the OVI UV 
doublet [OVI-UV]. These measurements, however, are at the 
limit of current instrumentation's capability. 

So far, the best WHIM detection has been provided by the 
analysis of the UV spectra of 31 AGN by Danforth & Shull 
(2005). They detected 40 OVI-UV absorption systems coun- 
terparts to Ly-a absorbers, whose cumulative number per unit 
redshift above a given equivalent width agrees with the latest 
theoretical predictions (Cen & Fang 2007). Most of the Ly-a 
lines associated to these absorbers, however, are too narrow 
to be produced by gas at temperatures larger than 10^ K. This 
means that either these lines are associated to the local Ly- 
a forest, and thus their contribution to the total baryon bud- 
gets has been already accounted for (Tripp et al. 2008), or 
that two different gas phases are found at the same location: 
a photo-ionized gas responsible for the Ly-a absorption and 
shock-heated WHIM responsible for the OVI lines. Assum- 
ing that all OVI-UV lines are indeed produced by a multi- 
phase gas, one obtains an upper limit to the mass density as- 
sociated to these absorbers which corresponds to 7-10 % of 
the missing baryonic mass (Danforth & Shull 2005 and Tripp 
et al. 2006). Alternatively, one could search for broad Ly- 
a lines in the UV spectra, i.e. lines with a width parameter 



6 > 40 km associated with gas hotter than 10^ K. The 
recent searches by Richter et al. (2006) and Stocke, Danforth 
& Shull (2008) showed that the few broad absorbers detected 
so far contribute to 2-10 % of the mass of the missing baryons. 

All these studies show that while some of the missing 
baryons have been revealed through OVI lines, the vast ma- 
jority of them is too hot to produce significant OVI opacity. 
Therefore the bulk of the WHIM can only be revealed and 
studied in the soft X-ray band. In emission, the diffuse signal 
detected in correspondence with overdensity of galaxies has 
been interpreted as the signature of the WHIM by Zappacosta 
et al. (2005) and by Mannucci et al. (2007). In absorption, 
among the several claimed WHIM detections at z > (e.g. 
Fang et al. 2002, Mathur et al. 2003), the most convincing 
one has been provided by Nicastro et al. (2005) who detected 
several WHIM absorption lines in the X-ray spectrum of the 
extraordinarily bright Blazar Mrk 421. Yet, the statistical sig- 
nificance of these detections has been questioned by Kaastra 
et al. (2006) through Monte Carlo simulations. Moreover, the 
absorption lines originally seen in the spectrum obtained with 
the Low Energy Transmission Grating Spectrometer on board 
of Chandra were not detected in a spectrum of the same object 
subsequently obtained through a longer exposure with the Re- 
flection Grating Spectrometer on XMM-Newton (Rasmussen 
et al. 2006). The very existence of this unsettled contro- 
versy demonstrates the importance of performing a dedicated 
experiment to detect the missing baryons in the X-ray band. 
Even more, going beyond the mere detection of the WHIM, 
the characterization of its thermal and chemical properties, 
will open an unexplored territory for the study of the interac- 
tion and co-evoultion of the diffuse phase of cosmic baryons 
and the stellar population. 

The possibility of detecting the WHIM in the absorption 
spectra of X-ray bright objects like AGN has been thoroughly 
discussed from a theoretical perspective in a series of stud- 
ies. While Chen et al. (2000) and Cen & Fang (2006), to 
mention a few examples, addressed the problem in a generic 
cosmological setup, other authors (e.g. Kxavtsov et al. 2002, 
Klypin et al. 2003, Viel et al. 2005) focused on the the gas 
distribution in the local Universe. In addition, other works 
have specialized the treatment to some specific satellites, like 
XMM-Newton (Kawahara et al. 2006) or the proposed DIOS 
(Yoshikawa et al. 2003, Yoshikawa et al. 2004), XEUS and 
Constellation-X (Viel et al. 2003) satellites. A comprehen- 
sive review of the current and future instrumentation for the 
study of the WHIM can be found in Paerels et al. (2008). In 
this work we focus on the possibility of revealing the missing 
baryons signature in the X-ray spectra of the GRB afterglows. 
The idea of using GRBs as bright X-ray beacons to study the 
physical properties of the intervening gas and trace its evo- 
lution out to very large redshifts has been first proposed by 
Fiore et al. (2000). Here we push this idea further and ex- 
ploit the recent SWIFT results on the rate of bright GRBs to 
estimate the probability to detect the WHIM and to character- 
ize its physical status in the context of the recently proposed 
EDGE and XENIA space missions. In addition, rather than 
relying on some specific theoretical prediction, we consider 
different WHIM models obtained from hydrodynamical sim- 
ulations to account for theoretical uncertainties. 

The characteristics of the detectors on board of EDGE 
and XENIA can be found in Piro et al. (2008), in 
Kouvelioutou et al. (2008) and at the web address 
http://projects.iasf-roma.inaf.it/EDGE Gamma-Ray Bursts 
are localized by a Wide Field Monitor, a hard X-ray detec- 
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tor covering about 1/4 of the sky. The position derived on 
board is used to command a satelHte autonomous slew, that 
enables the acquisition of the GRB location by the Wide Field 
Spectrometer, an X-ray telescope equipped with high spectral 
resolution TES microcalorimeters, in 60 s. 

Here we list for reference the capabilities of the two instru- 
ments that are more relevant to this work: 

• Wide Field Spectrometer [WFS]: Energy range 0.2 — 
2.2 keV. Energy resolution l\E = 3 eV (goal: 1 eV) 
at 0.5 keV (TES microcalorimeters). Field of view 
FOV 0.7° X 0.7°. Effective area A=1000 cm^ at 0.5 
keV. Angular resolution: /S.6 — 2 arcmin. 

• Wide Field Monitor [WFM] : Energy range 6-200 keV. 
Energy resolution l:^E/E = 3% at 100 keV. Field of 
view FOV > 2.5 sr. Effective area A=1000 cm^ at 
100 keV. Location accuracy: — 3 arcmin. 

This work is organized as follows. In Section 2 we describe 
our WHIM models and discuss the characteristics of the hy- 
drodynamical simulations that we have used to extract the 
mock X-ray absorption spectra. In Section 3 we analyze the 
statistical properties of the absorption lines in the mock X-ray 
spectra, characterize the physical properties of the absorbing 
material and investigate the statistical correlations among the 
absorption lines. In Section 4 we quantify the probability of 
detecting WHIM lines in the GRB afterglow spectra and es- 
timate the number of detections expected with a satellite like 
XENIA or EDGE. In Section 5 we investigate the possibility 
of studying the physical status of the WHIM and to measure 
its cosmological abundance. We discuss our main results and 
conclude in Section 6. 

2. MODELING THE WHIM 

The best WHIM models currently available are based on 
numerical hydrodynamical simulations and provide robust 
predictions for the cosmological WHIM abundance as a func- 
tion of redshift, the filamentary appearance of its large-scale 
distribution and its thermal state (Cen & Ostriker 1999a, 
1999b, Chen et al. 2001, Cen & Ostriker 2006 to pro- 
vide a few examples). The latter is mainly determined by 
hydrodynamical shocks resulting from the build-up of cos- 
mic structures at scales that have become nonlinear. Addi- 
tional, non-gravitational heating mechanisms that depend on 
the little known feedback processes following star formation 
in galaxies. However, they are relevant only in high density 
environments in which the role of radiative cooling must also 
be taken into account. In this work we do not rely on a sin- 
gle WHIM model that includes sophisticated recipes to treat 
feedback effects. Instead, in an attempt to bracket model un- 
certainties, we use two different hydrodynamical simulations 
to develop three WHIM models with different phenomenolog- 
ical treatments of the non-gravitational heating mechanisms, 
star formation and feedback effects. 

The first WHIM model (dubbed model M in Table 1) re- 
lies on the same simulations described in Pierleoni et al. 
(2008). This numerical experiment, of which we only con- 
sider the outputs at z < 1, was performed with the GADGET- 
2 Lagrangian code (Springel 2005) in a computational box 
of size Q{)h~^ comoving Mpc, loaded with 480"^ Dark Mat- 
ter [DM] and 480^ gas particles. The Plummer equivalent 
gravitational softening was set equal to 2.5 kpc h^^ in co- 
moving units for all particles. This choice of parameters rep- 
resents a good compromise between box size and resolution 



and allows to investigate the WHIM properties down to the 
Jeans scale. Radiative cooling and heating processes were 
followed according to Katz, Weinberg & Hernquist (1996). 
The simulation assumes the same X-ray and UV background 
at z = that Viel et al. (2003) used to determine the ion- 
ization state of the intergalactic gas. The X-ray background 
is not computed self consistently but is modeled a posteri- 
ori as Ix = lSiiE/Exy^-^^e^p{-E/Ex), with = 
1.75 X 10^26 jj^-i sr-i s-i and Ex = 40keV, in 

agreement with observations (Boldt 1987; Fabian & Barcons 
1992; Hellsten et al. 1998). The UV background is of the 
form luv = /&y(£;/13.6eV)-l■^ with I^y = 2.3 x lO^^a 
erg cm^^ Hz^^ sr^^ s^^ , which agrees within a factor of 2 
with the estimate of Shull et al. (1999) and shows no appre- 
ciable differences with the intrinsic Haardt-Madau like spec- 
trum (Haardt & Madau 1996) used in the simulation run. The 
cosmological parameters of this simulation are similar to the 
WMAP-1 year values (Spergel et al. 2003) and refer to a flat 
ACDM 'concordance' model with = 0.3, fib = 0.0463, 
Ho = 72 km s^^ Mpc"\ n = 0.95 and as = 0.85. In 
Fig. [U we show the gas distribution at 2 = of a slice of 
thickness 6 Mpc (comoving) extracted from the compu- 
tational box. 

The star formation criterion simply converts all gas parti- 
cles whose temperature falls below lO'"^ K and whose density 
contrast is larger than 1000 into (collisionless) star particles. 
No energy feedback from supernova explosions, metal ejec- 
tion and diffusion was considered. After determining the ther- 
mal state of each gas particle we assign the metal abundance 
in the post-processing phase according to a phenomenolog- 
ical metallicity-density relation, calibrated on hydrodynam- 
ical simulations in which metal enrichment is treated self- 
consistently. In model M we have assumed the determin- 
istic relation Z = min(0. 2, 0.02(1 + 5^^/^) of Aguirre et 
al. (2001), where Z is the gas metallicity in solar units and 
Sg = pg/ (pg) — 1 is the gas overdensity with respect to the 
mean, (pg). To model the WHIM properties we have used all 
available outputs of the hydro simulation with z < 1 using 
the stacking technique described in Roncarelli et al. (2006). 

The second simulation that we use has been carried out by 
Borgani et al. (2004) using the same GADGET-2 code. The 
box size is larger that in the previous case and consists of a 
comoving cube of 192 Mpc h^^ side containing 480'^H-480^ 
dark matter and gas particles. The Plummer-equivalent grav- 
itational softening is e=7.5 h^^ kpc at z = 0, fixed in phys- 
ical units between z — 2 and z — and fixed in comoving 
units at earlier times. The cosmological model is a flat ACDM 
with f^A = 0.7, = 0.04, Ho = 70 km s^^ Mpc"^ and 
fJs = 0.8. Differently from the previous case, this simulation 
includes a treatment of the main non-gravitational physical 
processes that influence the physics of the gas: /) the star for- 
mation mechanism, that is treated adopting a sub-resolution 
multiphase model for the interstellar medium (Springel & 
Hernquist 2003), //) the feedback from SNse including the 
effect of weak galactic outflows and the radiative cool- 
ing assuming zero metallicity and heating/cooling from the 
uniform, time-dependent, photo-ionising UV background by 
Haardt & Madau (1996). 

We have used all available outputs of the simulation with 
z < 1 to construct two more WHIM models that differ by 
their metallicity-density relation. The rational behind speci- 
fying the metal abundance in the post-processing phase rather 
than adopting the one computed self consistently in the sim- 
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X [ h ' Mpc ] 

Fig. 1. — The spatial distribution of the gas overdensity 5g in a slice of 
thickness 6 Mpc (comoving) at z = extracted from the simulations 
described in Pierleoni et al. (2008) and centered around the largest cluster 
in the simulation box (to be found in bottom-right part of the panel) 

ulation is that the simple metal diffusion mechanism in the 
SPH code reproduces the correct metallicity in high den- 
sity regions like the intracluster medium, but underestimates 
the metal abundance in low density environments like the 
WHIM. The first model obtained from the simulation of Bor- 
gani et al. (2004), dubbed i?l in Table 1, assumes the de- 
terministic metallicity-density relation of Croft et al. (2001): 
Z = min (0.3,0.005(1 + Sgf/'^)- Several hydrodynamical 
experiments (e.g. Cen and Ostriker 1999a, 2006), have re- 
vealed that the relation between Z and 5^ is characterized by 
a very large, non-Gaussian scatter around the mean, reflecting 
the fact that the thermal state and the metal content of the gas 
are not uniquely determined by the local density but, rather, 
by its dynamical, chemical and thermal history. In model B2 
we account for this stochasticity by enforcing in our hydrody- 
namical simulation the same scatter found by Cen & Ostriker 
(1999a, 1999b) in their numerical experiment. To do this we 
have measured the density-metallicity relation in the simu- 
lation outputs of Cen & Ostriker (1999a) and constructed a 
2D probability distribution in the Z ~ 6g plane by counting 
particles in each bin, at various redshifts. The resulting 2D 
grid provides a probability distribution that we implement in 
a Monte Carlo procedure to assign metallicity to the gas par- 
ticles in the Borgani et al. (2004) simulation with known 
density. 

The rationale behind considering three different WHIM 
models is to provide some indication of the theoretical un- 
certainties that reflect our little understanding of the effects of 
stellar feedback from star- and galaxy formation, chief among 
which is of course the metal enrichment. The latter is so 
important in determining the observational properties of the 
WHIM that model predictions can be obtained by specify- 
ing gas metallicity in the post-processing phase come close 
to those computed self-consistently in the simulations (Chen 
et al. 2001). In what follows we mainly focus on model 



TABLE 1 
WHIM Models 



Model 


Lbox 


Ndm - 


f Ngas 


€ 


Qb 


0-8 


Z-pg relation 


M 


60 


400^ - 


h400=* 


2.5 


0.046 


0.85 


ry 1/3 

Zoc pg' 


Bl 


192 


480^ - 


h4803 


7.5 


0.04 


0.8 


ry 1/2 

Zoc Pg' 


B2 


192 


480^ - 


h480=* 


7.5 


0.04 


0.8 


Z-pg + scatter 



Column 1: Model name. Column 2: Box size (in comoving Mpc h-^). 
Column 3: Number of DM and gas particles. Column 4: Plummet equivalent 
gravitational softening (in kpc h~^). Column 5: Baryon density parameter. 
Column 6: Mass vaiiance at 8 Mpc h~^. Column 7: Metallicity-Density 
relation adopted. 



B2 that we regard as more realistic. However, we also con- 
sider the two models Bl and M, which adopt ideal metal- 
licity models directly used in previous analyses that have ap- 
peared in the literature. We stress that these models do not 
bracket all unknows in the theoretical calculation since they 
adopt simple recipes for metal enrichment and stellar feed- 
back processes and ignore departures from equilibrium that, 
for example, have been considered in the more sophisticated 
WHIM model of Cen & Fang (2006). 

After having simulated the metals and their distribution in 
the box and after having specified the metagalactic UV and 
X-ray background, we have used the photoionization code 
CLOUDY (Ferland et al. 1998) to compute the ionization 
states of these metals under the hypothesis of hybrid coUi- 
sional and photoionization equilibrium. The resulting ioniza- 
tion fractions turns out to be very similar to that obtained by 
Chen et al. (2003) and shown in their Figs. 2 and 3 that im- 
plies that our OVII fraction in regions of moderate overdensity 
is different from that computed by Kawahara et al. (2006) 
using a modified SPEX code (see their Fig. 1). The density 
of each ion was obtained through: nj{x) — nH{x)XjYzQ'Zi, 
where Xj is the ion fraction determined by CLOUDY that 
depends on gas temperature, gas density and ionizing back- 
ground, nn is the density of hydrogen atoms, Z is the metal- 
licity of the element in solar units, and Yz^ is the solar abun- 
dance of the element (Asplund et al. 2005). 

The ions responsible for the most prominent X-ray absorp- 
tion lines are listed in Table 2together with the energy of the 
transition and the oscillator strength. These are the X-ray lines 
that we simulate in our mock absorption spectra. In addition, 
to compare model predictions with existing data we also con- 
sider the OVI UV doublet (at 1 1.95 eV, 12.01 eV), which has 
been observed in the absorption spectra taken by FUSE and 
HST (Tripp et al. 2000, Savage et al. 2002, Danforth & 
Shull 2005). 

To check the validity of the of hybrid ionization-colUsional 
equilibrium hypothesis, we have computed the typical recom- 
bination time in the WHIM tree = {nwHiM ctree)^^, with 
aree the recombination rate. For the two more relevant ions 
OVII and OVIII, in the temperature ranges considered here we 
findttrec ~ 2xl0~"cm^s"i (Mazzotta et al. 1998), where 
both the dielectronic and radiative recombinations are taken 
into account. This means that the Hubble time is larger than 
the recombination time for a density > 10~^ cm~^. Thus, for 
a large fraction of the WHIM, the approximation of photoion- 
ization equilibrium should be reasonable. Indeed, as it has 
been shown by Yoshikawa and Sasaki (2006), departures from 
equilibrium do not dramatically alter the equivalent widths 
[EW hereafter] of the OVIII and OVII absorption lines. 
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TABLE 2 
Simulated WHIM Ions. 



Ion 


Energy [eV] 


/ 


OVI 


11.95+12.01 


0.19 


OVI ls-2p 


563 


0.53 


OVII ls-2p 


574 


0.70 


OVII ls-3p 


666 


0.15 


OVIII ls-2p 


654 


0.42 


CV ls-2p 


308 


0.65 


CVI ls-2p 


367 


0.42 


NelX ls-2p 


922 


0.72 


FeXVII 2p-3d 


826 


3.00 


MgXI ls-2p 


1352 


0.74 



Column 1: Ion. Column 2: Wavelength [eV]. Column 3: Oscillator Strength. 
2.1. Mock absorption spectra 

Having determined the density of a given ion along the line 
of sight [LOS], the optical depth in redshift-space at velocity 

u (in km s^^) is 



Tj{u) = 



crpjc 
H{z) 



dyni{y)V u - y - v\{y),b{y) 



(1) 



where cto./ is the cross-section for the resonant absorption 
which depends on the wavelength. A/ and the oscillator 
strength // of the transition , y is the real-space coordinate 
(in km s^^), V is the standard Voigt profile normalized in 
real-space, h = {2kBT / mic^Y /'^ is the thermal width. Ig- 
noring the effect of peculiar velocities, i.e. setting ~ 0, the 
real space coordinate y and the redshift z are related through 
d\/\ = dy/c, where A = A/(l + z). Since WHIM absorbers 
have low column densities, the Voigt profile is well approxi- 
mated by a Gaussian: V = (-y7r6)^^ exp[— (w — The 
transmitted flux is simply T — exp(— r). 

For each WHIM model we have extracted 15 independent 
spectra from the light-cones obtained by stacking all available 
simulation outputs. Mock spectra from model M have a spec- 
tral resolution of 3 km and extend out to z = 0.5. Spectra 
from models B\ and B2 have a resolution of '-^ 9 km s^^ 
and extend out to z = 1.0. Figure |2] shows, for model Bl. 
and from top to bottom, the transmitted flux in OVII, OVIII, 
the density and the temperature along a generic LOS, out to 
z = 0.1. Strong absorption features are found in correspon- 
dence of density and temperature peaks. OVIII lines sample 
regions of higher density and temperature than OVII lines. 
The typical number of absorbers is small and they can be eas- 
ily identified in the mock spectra using a simple threshold al- 
gorithm. The positions of the absorption lines are identified 
with the minima of the transmitted flux and their EW are 
measured by integrating J-{u) over the redshift interval [z, z] 
centered around the minimum, in which T < Tthr- The few 
overlapping lines are separated in correspondence of the local 
maximum between the two minima. The threshold to identify 
absorption lines is Tthr = 0.95. We have checked that results 
do not change significantly if one lowers the threshold to 75 
% of the transmitted flux. 

3. LINE STATISTICS 

To check the reliability of our WHIM models we first com- 
pare our theoretical predictions with the only WHIM unam- 
biguous detections available so far In Figure [5] the differ- 
ent symbols connected by continuous lines show the cumula- 
tive distributions of OVI UV doublet EWs predicted by our 
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Fig. 2. — Simulated spectrum from model B\ for a generic LOS with two 
strong absorbers. From Top to Bottom: Transmitted flux for OVII, OVIII, 
temperature and density along the LOS as a function of redshift. 

models. Error bars represent the scatter among the different 
mock spectra. Filled triangles with error bars show the cu- 
mulative EW distribution of the OVI absorbers measured by 
Danforth & Shull (2005). Over most of the EW range, all of 
the considered WHIM models are consistent with data within 
the errors. The spread in model predictions reflects theoret- 
ical uncertainties. The larger number of OVI lines in model 
B2 results from the scatter in the density-metallicity relation 
that systematically increases the probability of an OVI absorp- 
tion line. Model AI predicts a cumulative distribution steeper 
than models Bl and B2. At large EWs the discrepancy re- 
flects the different amount of large scale power in the com- 
putational boxes of the parent simulations. At small EWs, 
instead, the difference originates from the different numerical 
resolution of the two simulations. In this work we focus on 
absorption lines with EW > 10 km s^^. Hereafter we use 
different units for the line EWg- Useful conversions between 
units are EW = 6.67(i;T^/100 km s-i)(A/20 A) mA = 
0.33{EW/100 km s-^){E/l keV) eV. 

The scatter in the model predictions for the OVI line distri- 
bution function is quite small and one may wonder whether 
this is also the case for the X-ray absorption lines. Let us 
consider the two strongest WHIM lines: OVII and OVIII. 
These ions probe regions of higher density and temperature 
than OVI. In such environments the predicted ion abundance 
is more sensitive to the different metallicity-density prescrip- 
tions adopted in the models and thus we expect a larger spread 
among model predictions. Figure |4] shows that this is indeed 
the case: the differences among the predicted cumulative dis- 
tributions for OVII (left panel) and OVIII (right panel) are 
larger than in the OVI case. Predictions of model B2, that 
account for the scatter in the metallictity, are consistent with 
those of the Chen et al. (2003) model, in which the metal- 
licity is drawn randomly from a lognormal distribution with 
(Z/Zq) = -1.0, (Tiogz = 0.4. Model 52 predictions ai-e 
also consistent with those of Cen & Fang (2006), that include 
the feedback effect of galactic super-winds and account for 
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Fig. 3. — Cumulative distribution function of OVI UV doublet equivalent 
widths. Model predictions vs. observations. Filled triangles with eiTor bars: 
Danforth & Shull (2005) measurements. Filled squares and dot-dashed green 
line: model M. Crosses and dashed red line: Model Bl. Filled dots and soUd 
blue line: Model B2. For each model the eiTor bars represent the scatter in 
the 15 independent mock spectra. 
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Fig. 4. — Cumulative distribution function of absorption lines EWb- 
Model predictions. Left: cumulative distribution function of OVII. Right: cu- 
mulative distribution function of OVIII. Filled squares and dot-dashed green 
line: model M. Crosses and dashed red line: Model Bl. Filled dots and 
solid black line: Model B2. Error bars in the model represent the scatter 
among the 15 independent mock spectra. 

departures from local thermodynamic equilibrium. We note 
that all our predictions, including those from models M and 
Bl, satisfy the current upper limits on the OVII and OVIII 
line observations. 
In Figure |5] we show, for model B2 only, the predicted cu- 
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Fig. 5. — Cumulative distribution function of absorption lines equivalent 
widths. Predictions of Model B2. Left: cumulative distribution function of 
OVI X-ray (blue continuous) OVII ls-3p (red long-dashed), CV (green dot- 
ted) and CVI (black short-dashed). Right: cumulative distribution function 
of FeXVII (blue solid), MgXI (red long-dashed), and NelX (black dotted). 

mulative distribution functions for the other ions listed in Ta- 
ble 2. As already pointed out by Chen et al. (2003) CVI and 
CV are the strongest absorption lines after oxygen. However, 
their transition energies of 367 and 308 eV put them in a range 
of poor instrument sensitivity (see Fig. [Tol ) and restrict their 
observability to low redshift sources. NelX produces the next 
strongest absorption together with the OVII ls-3p line. The 
latter, however, is observed in association with the stronger 
OVII ls-2p line and, as we will discuss in the next Section, 
improves the statistical significance of the WHIM detection, 
but offers no further information on the absorber's physical 
state. On the contrary, the OVI line in the X-ray band could, 
together with OVII and OVIII, be used to determine unam- 
biguously the physics of the gas. Unfortunately the X-ray 
OVI line is typically ^ 4 times weaker than the OVII line and 
will be difficult to detect unless instrumental sensitivity will 
increase dramatically. For this ion UV spectroscopy appears 
as the only effective option currently available. 

3.1. Physical properties of mock absorbers 

Thanks to our WHIM models we can substantiate our pre- 
vious statement that observed OVI lines have probed a few 
per cent of the missing baryons that reside in the relatively 
cool regions. In Fig.|6]we show the phase-space diagram of 
the cosmic baryons at z in the hydrodynamical simula- 
tion of Borgani et al. (2004). Small-size points represent gas 
particles in the simulation. Some regions of this phase-space 
diagram have been already probed by observations, like the X- 
ray emitting gas in the central regions of virialized objects (the 
green points with Tg > 10^-^ K and pg/(pg) > lO^-^), the dif- 
fuse, cold gas responsible for the local Ly a forest (blue points 
with Tg < lO'^ K and Pg/{pg) < 10^) and the cold, dense gas 
observed in association to star forming regions within galax- 
ies (yellow points with Tg < 10^ K and Pg/(pg) > 10^). The 
majority of the remaining baryons are located in regions of 
the phase space that are currently not probed by observations. 
Some of them are expected to be found in the outskirts of the 
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galaxy clusters all the way out to the virial radius (red points 
with Tg < lO^-S K and lO^-^ < Pg/(pg) < lO^-^) but the 
large majority, which constitute the bulk of the WHIM, is to 
be found in lower density environments 

The larger symbols characterize the gas that can be detected 
through OVI (filled white squares), OVII (filled magenta dots) 
or OVIII (filled cyan triangles) assuming model B2. Only 
lines with EW > 10 km are considered in the plot. Dif- 
ferent ions preferentially probe different regions of the phase 
space with a non negligible overlap, spanning a large range of 
WHIM physical conditions. As anticipated, OVI is a signpost 
for the cooler regions of the WHIM while, on the opposite, the 
OVIII lines typically flag hotter environments. The OVII ab- 
sorbers span a very large range in temperature, although those 
with Tg < 10^ K are typically associated with weak lines. 

The position of the OVI, OVII and OVIII symbols in the 
phase space is not determined precisely since the density and 
the temperature of the gas associated to the absorption line are 
not uniquely defined. In this work we estimate temperature 
and density of the absorbing material as in Chen et al. (2003). 
First we define the temperature and density in each spectral 
bin as 

Pg(") = -TT / ^Pgiy)d'y; Tg{u) = [ ^Tg{y)dy 

t(u) J dy t{u) J dy 

(2) 

and then we obtain the corresponding quantities for the ab- 
sorbers by averaging over the bins in the line redshift interval 
[z , 1] . This is the definition that we have used to plot the sym- 
bols in Fig.|6] 

We have also considered an alternative estimator in which 
temperature and density are weighted by the optical depth of 
the line: 

^ _ Il'r{y)pgiy)dy ^ ^ J^T{y)Tg{y)dy 

^ Sl'r{y)dy ' ° JlT(y)dy 

Using Eq.[3]rather Eq.|2]does not introduce systematic differ- 
ences, just a random scatter of ^ 15% in both temperature 
and density. 

Another possible source of ambiguity in estimating the tem- 
perature ad the density of the gas is represented by those ab- 
sorption systems characterized by two or more ions, each one 
with its own opacity. In this case, the redshift interval [z, 1] 
depends on the ion considered and so will be the estimated 
temperature and density of the absorbing gas. To assess the 
importance of this effect we have considered all mock ab- 
sorption systems featuring any two lines among OVI, OVII 
or OVIII and found ^ 20 % scatter in temperature and den- 
sity computed in the redshift ranges of the different lines. 

3.2. Correlation between absorbers 

Figure |2] shows the predicted correlation between the equiv- 
alent widths of the oxygen lines associated to the same ab- 
sorber in the mock spectra of model B2. We only consider 
absorbers in which all lines are detected above threshold. To 
avoid ambiguity we measure the line EW in the velocity 
range used to identify the OVII line. In the plots crosses indi- 
cate lines that are not saturated (i.e. for which the transmitted 
flux is above zero). Green triangles (blue squares) indicate 
absorbers in which the ion line on the X (Y) axis is saturated. 
Red dots indicate that both ion lines are saturated. The bottom 
panel shows the correlation between the equivalent widths of 
the OVI UV-doublet and the OVI X-ray line. For optically 
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Fig. 6. — Phase space diagram of the cosmic baryons aX z < 0.2 from 
the hydrodynamical simulation of Borgani et al. (2004). Blue dots with 
T < 10^ and p/(p) < 10'^: gas particles detected in the local Ly a for- 
est. Yellow dots with Tg < 10^ K and pg/(pg) > 10'^ gas particles 
associated with star forming regions. Green dots with Tg > lO^"'^ and 
Pe/{Pg) ;^ 10^'^ gas particles in the central part of virialized structures 
that are potentially observable with current X-ray satellites. Red dots with 
Tg < lO^-'' and lO^'^ < Pg/(pg> < 10^'^ gas particles in the outskirts of 
virialized objects cun'ently not observed. Black dots: currently unobserved 
intergalactic medium including the warm-hot phase. Filled white squares gas 
detectable through OVI absorption lines. Filled magenta dots gas detectable 
through OVII absorption lines. Filled cyan triangles gas detectable through 
OVIII absorption lines. Only absorption lines with EW > 10 km s^^ have 
been considered. 

thin lines there is a simple, deterministic relation between the 
EWs of the two lines. The scatter in the plot arises from hav- 
ing measured the EW of the lines in the velocity range ap- 
propriate for the OVII line. Equivalent widths of saturated 
UV-absorbers deviate from the expected correlation, causing 
the flattening at large EWs- There is a small range of tem- 
peratures in which both OVII and OVI have high fractional 
abundance. In this interval both species trace the underly- 
ing baryon density, causing the tight correlation seen in the 
third panel of Fig|7] Strong OVII lines, however, have a large 
range of OVI EWs- Part of this spread is due to line satura- 
tion and part of it is due to the fact that in warm absorbers 
the OVI fraction is very low. As pointed out by Chen et al. 
(2003) who performed the same correlation analysis, this fact 
suggests that, at least for the stronger lines, follow-up X-ray 
observations of known OVI absorbers are a promising way 
to find OVII absorption lines while the opposite strategy is 
less efficient. Similar considerations apply to OVII and OVIII 
lines. The scatter in the relation, however, is smaller than in 
the previous case since these ions have large fractional abun- 
dance over a large range of temperatures, a fact that is clearly 
illustrated in the phase-space diagram of Fig|6] Finally, there 
is little or no correlation between OVI and OVIII except for 
strong lines. In this case the large scatter reflects the fact that, 
for a fixed gas density, the temperature of these two species 
can span a three order of magnitudes range. As a result, OVI 
lines are not able to efficiently flag the presence of OVIII ab- 
sorbers and vice versa. 
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Fig. 7. — Equivalent width correlation for oxygen ions. From top to bot- 
tom: OVIII vs. OVI-UV, OVIII vs. OVII, OVII vs. OVI-UV, OVI-UV vs. 
OVI-X. Black crosses indicate unsaturated lines. Filled green triangles: sat- 
urated line for the absorber on the X-axis. Filled blue squares: saturated line 
for the absorber on the Y-axis. Filled red circles: both ion lines saturated. 

4. DETECTING THE WHIM THROUGH ABSORPTION 

In this Section we explore the possibility of detecting 
WHIM in the absorption spectra of X-ray bright objects. We 
will focus on GRBs, showing that these transient sources offer 
several advantages over bright AGN like quasars and Blazars. 
Then we use our WHIM models to estimate the number of de- 
tections that one may expect from a XENIA/EDGE-like satel- 
lite mission. 

4. 1 . GRBs as background sources 

Unambiguous WHIM detections require X-ray background 
sources that are bright, distant and common enough to provide 
a good statistical sample of WHIM lines. The latter require- 
ment excludes rare objects, like the Blazars during their ac- 
tive phase. Bright AGN are quite common but nearby, hence 
probing relatively short lines of sight. Cross-correlating the 
VERONCAT catalog of Quasars and AGN (Veron-Cetty and 
Veron 2006) with the ROSAT All-sky Survey Bright Cata- 
logs (IRXS, Voges et al. 1999), Conciatore et al. (2008) 
found only 16 objects (mostly BLLac) at z > 0.3 with a 200 
ks integrated flux (fluence) in the [0.3-2] keV range above 
10^^ erg cm^^, bright enough to detect WHIM lines The 
mean redshift of these objects is r^ 0.5, to be compared with 
the redshift distribution of the Warm-Hot gas (see e.g. Fig. 
1 of Cen & Ostriker 2006) that, at z ~ 2 still accounts for 
~ 20% of the barionic mass in the universe. GRBs look more 
promising in all respects. First of all their X-ray afterglows 
provide a fluence in the [0.3-2] keV similar to that of bright 
AGN, collected however in only 50 ks rather than 200 ks. The 



only disadvantage is represented by their transient nature and 
by the fact that, since they emit a considerable fraction of the 
soft X-ray photons soon after the prompt, they require an in- 
strument capable of fast re-pointing. 

In this work we will focus on GRBs as background sources 
since the Wide Field Monitor proposed for the EDGE/XENIA 
missions is especially designed to trigger on and localize 
GRBs and X-ray Flashes and the satellite can re-point the 
source in 60 s. Spectra would be taken by the Wide Field 
Spectrometer in the 0.2 - 2.2 keV band by integrating the 
X-ray flux between 60 s and 50 ks. To convert counts to 
fluxes, or fluxes from one band to another, and to correct 
for absorption in our Galaxy and in the host galaxy we as- 
sume a typical power-law GRB spectrum with photon index 
r = 2, absorbed by Nh (Galaxy) = 2 x lO^o cm^^ and by 
NH(Sourcc) = 3.2 x IG^i cm^^. 

Given that EDGE/XENIA is expected to have orbit and slew- 
ing capabilities similar to those of Swift, a reliable way to es- 
timate the expected afterglow fluence distribution in a given 
energy band is to use the total number of photons actually de- 
tected for each GRB by Swift/XRT in that energy band and 
converting counts to physical units by adopting the best fit 
spectral model of each GRB. In this way, all the effects that 
may affect the fluence measured by the X-ray telescope (e.g.: 
different times to go on target. Earth occultation statistics and 
South Atlantic Anomaly passage for a satellite in Low Earth 
Orbit, combined with the different shape and duration of both 
X-ray prompt and afterglow emissions. X-ray flares, long du- 
ration prompt emission, possibility of triggering on a precur- 
sor and thus being on target during the main phase of prompt 
emission, etc.) are automatically taken into account. 

In the light of these considerations, we have considered, for 
a sample of 187 GRBs, the total number of counts measured 
directly by Swift/XRT in the [0.3-1] keV band during foUow- 
up observations. We stress that the fluence in the narrow [0.3- 
1] keV energy band allows to trace well the monochromatic 
0.5 keV fluence which is more relevant to assess the possi- 
bility of detecting WHIM lines. In fact, the errors caused by 
the possible scatter of the GRB population around the typical 
spectrum, which is used to convert the observed counts in a 
0.5 keV fluence (see column 1 in Table 3), are minimized by 
using the narrower [0.3-1] keV energy band. 

Based on this approach, we estimated that the afterglows 
with [0.3-1 keV] total counts corresponding to a monochro- 
matic fluence 0.5 keV higher than 8.8, 22 and 43x 10^^ erg 
cm^^ keV^^ are, respectively, 9%, 3% and 1% of the to- 
tal. These fluence thresholds are listed for reference in the 
first column of Table 3. In the second column we also list 
the corresponding value of the [0.3-10] keV unabsorbed (i.e., 
corrected for the absorption of the Galactic and the intrinsic 
hydrogen-equivalent column density) fluence, derived from 
the [0.3-1] keV counts by assuming the typical GRB after- 
glow spectrum specified above. Taking into account the Swift 
detection rate of r^lOO GRBs per year and the larger solid an- 
gle of the WFM on board of EDGE/XENIA with respect to the 
SWIFT/BAT monitor (a factor ^2 for bursts with prompt flu- 
ence of about 1 X 10^® erg cm^-^, and an additional factor ^^^3 
for larger values of the fluence, see Piro et al. (2008) for 
details), one expects to observe about 18, 8 and 3 such events 
per year, respectively. 

These numbers represent a conservative estimate based on 
current Swift capabilities. They can be increased in two ways. 
First of all, reducing the slewing time from the average 120 
s required by Swift to 60 s would increase the count fluence 
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Fig. 8.— X-ray [0.3-10] keV fluence distribution for 125 GRBs observed 
by Swift with observation start time < 1000 s. The three histograms assume 
different start times and uninterrupted integration intervals. Red curve: 60 s 
- 60 ks. Blue curve: 120 s - 60 ks. Black curve limiting the shaded area: 180 
s - 60 ks. 



by a factor ~2. The ampHtude of this effect is quantified 
by the change in the distribution of the unabsorbed fluence 
in the [0.3-10] keV shown in Fig. [S] where different obser- 
vation start times from the burst onset are assumed. These 
fluences have been derived from a sample of 125 Swift/XRT 
light curves, assuming uninterrupted observations (the ob- 
served flux in Swift, due to the low earth orbit, is basically 
half of that of an uninterrupted observation), and converting 
counts to physical units by adopting the best fit spectral model 
of each GRB. The blue histogram shows the case of a start 
time from the burst onset of 120 s, equal to the average Swift 
case. Using a satellite like EDGE/XENIA, for which the slew- 
ing time could be reduced to 60 s, would shift the distribution 
toward significantly larger values of the fluence, as shown by 
the red histogram. In that histogram there are 3 events exceed- 
ing 10^'' erg cm^-^. Two of them are due to extrapolation of 
the initial very steep decay. The third one is a really bright 
but very peculiar event (GRB 060218) with a very long (> 35 
min) prompt emission which was seen by Swift/XRT too. 

The second improvement is represented by the possibility 
of detecting and following up of the softest GRB, i.e. the 
X-Ray Rich (XRR) GRBs and the X-Ray Flashes (XRF) 
whose existence has been revealed by BeppoSAX (Heise et 
al. 2001). Their abundance turned out to be comparable to 
that of classical GRBs. In addition, the observations by Bep- 
poSAX (D'alessio et al. 2006) recently complemented by 
Swift XRT (Gendre et al. 2007) showed that, despite of a 
less energetic prompt gamma-ray emission, the X-ray after- 
glow flux of XRFs and XRRs is similar to that of GRBs. As 
a consequence, detecting these populations of objects would 
significantly increase (almost by a factor of two) the number 
of bright afterglows that can be used to study the WHIM. In- 
deed XRFs constitute about 30% of the HETEII sample, while 
the combined XRFh-XRR sample amounts to 70% of the to- 
tal (Sakamoto et al. 2005). For reference, XRFs constitute 
only about 10% of the Swift sample (Gendre 2007). This dif- 
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Fig. 9. — Cumulative distribution of the GRB X-ray afterglows detected 
per year by a satellite like EDGE/XENIA as a function the unabsorbed 
fluence in the [0.3-1] keV band. Vertical lines are drawn at the three different 
fluence thresholds of 2, 5 and 10 X lO^'' erg cm^^ in the [0.3-10] keV band. 
These have been obtained assuming the typical GRB afterglow spectrum and 
column density described in the text. 

ference is due to the higher lower energy threshold (15 keV) 
of the BAT trigger instrument, compared to HETEII WXM (2 
keV) or to BeppoSAX WFC. Adopting a low energy thresh- 
old of 8-10 keV, which represents the baseline for the monitor 
planned for EDGE/XENIA, is expected to increase by ^30 % 
the number of detections. With a low energy threshold of 5 
keV the increase would be of ~ 70%. 

The cumulative distribution of GRBs that EDGE/XENIA 
is expected to observe in one year as a function of the 
counts measured in the [0.3-1] keV band is shown in 
Fig. |9] (see the EDGE science goal document available at 
http://projects.iasf-roma.inaf.it/ for further details). Assum- 
ing a typical GRB spectrum, we draw three vertical lines at 
the same reference 0.5 keV monochromatic fluences of 8.8, 
22 and 43 x 10^^ erg cm^^ keV^^ considered above, that cor- 
respond to the [0.3-10] keV unabsorbed fluence values listed 
in column 2 of Table 3. The number of expected detections 
above these reference fluence values is 48, 18 and 8 respec- 
tively. 

In column 3 of Table 3, we report the number of GRBs 
that EDGE/XENIA is expected to detect above a given fluence 
threshold, computed by averaging between the pessimistic 
case, obtained by simply scaling the Swift results, and the 
optimistic case, which accounts for the optimization in the 
slewing capabilities and the possibility of detecting XRFs and 
XRRs in addition to classical GRBs. 

As it can be seen, EDGE/XENIA is expected to detect in 
3 years ~ 100 GRBs above the fluence threshold indicated 
in Table 3 compared to the 16 distant AGN potentially useful 
for WHIM detection, if observed for 200-ks each. In fact, as 
we have stressed already, the situation is even more favorable 
for the WHIM detection since the average redshift of GRBs 
is significantly larger that AGN and the probability of line 
detections is expected to increase accordingly. 

4.2. Probability of OVII detection 

In this Section we quantify the detectability of WHIM ab- 
sorbers by applying the treatment and formalism of Sarazin 
(1989) to transient sources. We first focus on the problem of 
detecting the OVll line, which, as stressed before, is usually 
the sti'ongest absorption feature for the warm-hot intergalactic 
gas. 
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TABLE 3 
GRBS AND WHIM LINE DETECTION 



Fo.5 


F03-10 


^GRBs 


EWovi, 


^OVII 


EWoviii 


^OVII+OVIII 


8.8 10-'^ 


2 


33 


0.28 


58 


0.19 


37 


2.2 10-*^ 


5 10-6 


13 


0.18 


47 


0.12 


29 


4.3 10-*^ 


1 10-5 


6 


0.12 


32 


0.08 


19 



Column 1: Fluence at 0.5 keV [erg cm-^ keV-^]. Column 2: [0.3-10] 
keV Unabsorbed fluence [erg cm— ^] collected during the ORB follow-up 
observations. Column 3: Number of GRBs per year above fluence in a field 
of view of 2.5 sr. Column 4: Minimum EW (in eV) for a 4.5ct OVII line 
detection in the B2 model. Column 5: Number of 4.5(T OVII line detections 
per year in the B2 model. Column 6: Minimum EW (in eV) for a 3(7 OVIII 
line detection in the B2 model. Column 7: Number of simultaneous OVII 
and OVIII line detections per year with significance > 5 cr in the B2 model. 
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Fig. 10. — Mock X-ray absorption spectrum of a typical GRB afterglow 
at ^ = 2.6 convolved with the EDGE/XENIA response matrix. Dotted 
vertical lines mark the most prominent absorption lines produced by a WHIM 
filament at 2 = 0.069. A second absorption system is also present along 
the LOS at z = 0.298. Its absorption lines are flagged with the vertical 
dashed lines. The ions responsible for the absorption features are indicated in 
correspondence of each line. The fluence at 0.5 keV of the afterglow is also 
shown. 



In Fig.[TO]we show one realistic mock absorption spectrum 
of a GRB X-ray afterglow convolved with the EDGE/XENIA 
WES response matrix. Two WHIM absorbers at two different red- 
shifts (indicated in the plot) are present. In both cases the OVII 
absorption line stands out prominently along with the OVIII line. 
The nearest absorber at z = 0.069 also shows a weak CVI line. 
The absorber at z — 0.298 instead, is characterized by consider- 
ably stronger absorption features including the EeXVII, NelX and 
the OVII ls-3p lines. 

In order to detect an absorption line with statistical significance 
a, the detector must collect, within the instrument energy resolution, 

AE, a number of continuum photons Nphotons^ i'w)'^ ("^W) ^ ■ 
The corresponding minimum equivalent width EWmin for a given 



line fluence F is 

EWn,in 

100 km/s 
E 



1.3 X 10 ^ergcm 



IkeV 



AE 
10 eV 



-2koV"\ 
A 

\ 10000 cm2 



X 

-0.5 



(4) 



where A is the effective collecting area. 

In column 4 of Table 3 we list the minimum equivalent width re- 
quired to detect the OVII line with a significance of 4.5 a assuming 
an energy resolution AE — 3 eW and a collecting area A = 1000 
cm^ (i.e. matching the baseline configuration of the Wide Eield spec- 
trometer of EDGE/XENIA) in the X-ray spectra of a GRB with a 0.5 
keV fluence indicated in column 2 of the same Table. These esti- 
mates have been obtained by setting E = 0.574/(1 + 0.32) = 0.435 
keV in eq. |4]i.e. by assuming that all absorbers are at 2 = 0.32 
which is the average redshift of mock absorbers with EW > 0.08 
eV in the B2 model The reason for considering 4.5 a significance 
is that it guarantees that when two lines (typically OVII and OVIII) 
can be detected, the underlying absorber is detected with 5 a signif- 
icance, as will become clear in the next Section. The corresponding 
number of OVII line detections with significance > 4.5 a is listed 
in column 5. It was obtained by counting the number of OVII ab- 
sorbers with EW > EWn^in (in column 4) in the mock spectra of 
model B2 along a number of lines-of-sight equal to that of the ex- 
pected GRBs listed in column 3. We recall that the average number 
of OVII lines above EW is shown in the left panel of Eigure|4l in 
which we can appreciate the differences among the various WHIM 
models' predictions. Using models M or Bl instead of model B2 
which we regard as better physically motivated, significantly reduces 
the number of expected detections. The exact magnitude of the ef- 
fect depends on the EW of the line and on the GRB fluence. We 
estimate that, even with the most conservative model, a satellite like 
EDGE/XENIA should be able to detect about 25 OVII lines per year. 

We point out that, in all model explored, the number of expected 
OVII detections listed in Table 3 is biased low since in our mock 
spectra we have only considered OVII lines at z < 1 (z < 0.5 in 
model AI) while, as we have already pointed out, significant amount 
of WHIM gas are expected out to significantly higher redshifts. In 
addition, all these estimates should be boosted up by ~ 40 % if the 
energy resolution of the spectrograph could improve to 1 eV, which 
constitutes the goal for the WES. 

These estimates for a single OVII line detection, however, must 
be regarded as hypothetical since with no a-priori knowledge of the 
absorber's redshift there is no efficient way of preventing contami- 
nation by spurious absorption features or confusion with absorption 
lines form other ions. Of the two effects, confusion should be re- 
garded as the minor one given the paucity of expected non-Galactic 
absorption lines in the energy range relevant for the WHIM. An ef- 
fective procedure to minimize confusion is suggested by the analysis 
of our mock spectra. The idea is to identify the strongest absorp- 
tion feature with the OVII line (which is often the case), look for 
all other possible ion lines at the same redshift as the OVII line then 
move to the next strongest line and so on. The second effect, con- 
tamination by Poisson noise, should be dominant. However, its im- 
portance decreases with the EW of the OVII line. To quantify the 
amount of contamination we have produced a number of mock GRB 
spectra consisting in Monte Carlo realizations of power-law spectra 
with a specified Nh absorption, convolved with the EDGE/XENIA 
response matrix, all of them with a line fluence E0.5 in the range 
[0.9 — 4.3] X 10^*' ergcm-^ keV-^. These mock spectra have been 
"observed" to compute the probability of fake detections resulting 
from Poisson fluctuations as a function of the EW of the spurious 
line. By searching for OVII lines in the energy range [0.287—0.574] 
keV, corresponding to a redshift range z = [0 — 1], we computed 
the minimum EW required to guarantee a contamination level be- 
low 1%. The resulting EWs are the same as those listed in column 
4 of Table 3, i.e. they correspond to the minimum EWs required to 
guarantee a 4.5 u detection of the OVII line. In other words, a 4.5 
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a detection of a putative OVII line guarantees that contamination by 
Poisson noise is below 1 %. 

This exercise illustrates that the simple search for a single OVII 
line in the GRB spectra observed by EDGE/XENIA would provide 
25-80 line detections per year, where the lower bound is obtained 
by considering the most pessimistic among our WHIM models and 
the upper bound assumes a goal energy resolution of leV. We stress 
again that out of this sample, we expect at most 1 % to be spurious. 

Finally, we note that the smallest minimum EW in Table 3 for 
a joint OVII and OVIII lines detection, which we investigate in the 
next Section, is 0.08 eV. These values require a control of the sys- 
tematic effects, in particular of the energy width of instrumental bin 
and the level of the continuum, at ~ 2 % or better, corresponding to a 
systematics of 0.06 eV for 3 eV resolution (~ 8 % for a goal energy 
resolution of leV). 

4.3. Probability of multiple line detections 

X-ray spectroscopy allows in principle to characterize the physical 
state of the WHIM and measure its metal content. Unfortunately, the 
current energy resolution of the microcalorimeters in the soft X-ray 
band does not allow to resolve the WHIM lines and estimate the 
gas temperature directly from the line width. Indeed, typical line 
widths of detectable WHIM absorbers in our mock spectra are in the 
range 0.1-0.2 eV, well below the goal energy resolution of 1 eV of 
TES microcalorimeters. Nevertheless, as we shall discuss in the next 
Section, we can still use unresolved lines to constrain the physical 
state of the absorbing gas, provided that we can detect the absorption 
line of some other ion in addition to OVII. 

To asses the possibility of multiple line detection we first consider 
the OVIII line which often constitutes the second strongest absorp- 
tion feature in the mock spectra. Having already identified the red- 
shift of the absorber from the OVII line makes the detection of the 
OVIII line comparatively easier since one can now be happy with a 
3a rather than a 4.5(t detection. The minimum EW for such detec- 
tion can be obtained from Eq.|4]and is listed in column 6 of Table 3 
as a function of the GRB fluence. Like in the previous case of single 
OVII line detection, we have only considered OVIII absorption sys- 
tems with z < 1 and evaluated the minimum equivalent width of the 
line atz = 0.32, i.e. using E = 0.495 keV in Eq.H 

We find that ~ 50 % of all 3(7 OVIII lines are associated to a 
4.5(7 OVII line, i.e. the cumulative distribution function of OVIII 
and OVII line pairs lies a factor of two below that of the OVIII line 
alone, plotted in the right panel of Figure |4] The average ratio be- 
tween EW°y™ and BW'^in" is ~ 0.7, that is to detect the OVIII 
line associated to a 4.5 (7-detected OVII line the following condi- 
tion must be satisfied iSM^oviii > 0.7 -BWovii. Moreover, using 
Monte Carlo simulations we have verified that a joint OVII 4.5 a 
and OVIII 3 a detection corresponds to an overall ~ 5 a detection 
of the underlying absorption system, further reducing the probability 
of spurious lines contamination. 

Since the bulk of the WHIM is not in collisional ionization equi- 
librium, the EW ratio of two lines depends on both the temperature 
and the density of the absorbing gas. Therefore, the simultaneous 
observation of the OVII and OVIII lines alone cannot constrain the 
physical state of the WHIM. Instead, as we will show in the next Sec- 
tion, a third line, preferably of the same element, must be detected 
as well. To estimate the number of absorbers featuring at least three 
lines we looked for all mock absorption systems in which an addi- 
tional 3 (7 line is detected in association with a 4.5 a OVII line and 
3 (7 OVIII line. Table 4 shows the expected number of such simulta- 
neous detections per year as a function of the GRB fluence. We did 
not include OVII ls-3p in the Table since its detection would only 
increase the statistical significance of the absorption system but give 
no insight on its physical state (the EW ratio of the OVII ls-3p and 
OVII ls-2p lines only depends on their relative oscillator strength). 

Detecting 3 lines of the same elements, like OVI, OVII and OVIII, 
would allow to probe the physics of the WHIM directly. Unfortu- 
nately, as shown in Table 4, the X-ray OVI line is very weak and 
the number of expected multiple oxygen lines detection per year is 



TABLE 4 

Multiple WHIM lines detections involving OVII 1 s-2p. OVIII 

AND A THIRD ION 



^0.5 


^OVI-X 


#cv 


#CVI 


^NelX 


^FeXVII ^MgXI 


8.8 10-'^ 


4 


13 


19 


23 


7 6 


2.2 10"'' 


3 


12 


18 


23 


6 5 


4.3 10-'^ 


3 


7 


15 


15 


6 5 



Column 1: GRB fluence at 0.5 keV [erg cm"^ keV^^]. Column 2: Detec- 
tions per year involving OVI ls-2p. Column 3: Detections per year involving 
CV. Column 4: Detections per year involving CVI. Column 5: Detections 
per year involving NelX. Column 6: Detections per year involving FeXVII. 
Column 7: Detections per year involving MgXI. 



very small. For this purpose one should use the OVI UV line whose 
detection in the GRB spectra would require simultaneous UV and 
X-ray spectroscopy. 

Table 4 shows that the best chances to observe multiple line ab- 
sorption systems would be through the CVI or NelX lines. In these 
cases, however, to assess the physical state of the WHIM, the relative 
abundance of the elements needs to be specified a priori Sticking to 
these ions we have found very few cases of multiple detections which 
does not include both the OVII and the OVIII line. Also, a significant 
fraction of multiple line detections consists of more than three lines, 
i.e. the total number of multiple detections expected per year is less 
than the sum of all entries in the columns of Table 4. 

5. PROBING THE PHYSICAL STATE OF THE WHIM AND ITS 
COSMOLOGICAL ABUNDANCE 

The analysis of X-ray absorption spectra can provide much more 
than a mere WHIM detection. In this Section we address the problem 
of determining the temperature and density of the warm hot gas and 
its cosmological density. To explore the different WHIM models 
available, all results shown in this Section were first obtained for 
model Bl and then checked their validity for model B2. Because 
of the paucity of multiple-line absorption systems, we did not extend 
the analysis to the mock spectra of model AI. 

5.1. Measuring the density and the temperature of the WHIM 

As already mentioned, the quantitative analysis of the absorption 
spectra is complicated by the fact that the X-ray absorption lines 
cannot be resolved by currently proposed micro-calorimeters. On 
the other hand the weak WHIM lines are expected to be optically 
thin so that a simple linear relation holds between the line equivalent 
width and the ion column density: 

Wp(kms"') = \(X')}x'Nx^, (5) 

meC 

where X' represents the ion of the element X, Nxi = / ?ix> {y)dy 
its column density, A(Jf ) is the wavelength of the transition and 
fxi is the oscillator strength. From now on we use the subscripts p 
to indicate predicted values and m to indicate the values measured 
from the spectra. In Eq.[5]i?Wp is expressed in km s^^, the same 
units that we will use throughout this Section. 

Eq. [5] allows to estimate the column density of the ion from the 
EW of the absorption line. However, with no a-priori information 
on the size of the absorbing region, one cannot use the measured 
EW to estimate the actual density of underlying gas and assess its 
cosmological density. However, some insights on the physical state 
of the absorber can be obtained from the EW ratios of at least three 
different ion lines. In the optically thin limit the EW ratio between 
any two lines is: 

EWAX^^Wnfx^XlAx 
EWp{Yi) \(Yi) fyj Ay' 

where Xi and Yj represent the two ion fractions while the ratio be- 
tween Ax and Ay represents the relative abundance of the element 
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X compared to Y . Eq.|6]has been derived from Eq.[5]assuming that 
the density of the ions is constant across the absorber and can be used 
to infer the density and temperature of the gas upon which the ion- 
ization fractions Xi and Yj depend. Here we apply Eq.[6]to estimate 
the density and temperature of the absorbing systems in our mock 
spectra, hence using an approach similar to that of Nicastro et al. 
(2005) to analyze the physical state of the two putative WHIM ab- 
sorbers along the LOS to Mrk 421. For this purpose we use the same 
grids of hybrid coUisional ionization plus photoionization models 
that we have used to draw the mock spectra from the hydrodynam- 
ical simulations. These models are used to predict the ratio X^ /Y^ 
which we substitute in Eq.|6]to obtain EWp{X^) / EWp{Y^) as a 
function of gas density, temperature and ionization background for 
a given relative metallicity Ax /Ay, assumed a priori. In a sec- 
ond step, we compare the predicted value of EW-p{X^) / EWviY-') 
with that measured from the mock spectra, EWm (X' ) / EWm (Y^), 
to estimate the density and the temperature of the absorbing gas. Un- 
like Nicastro et al. (2005) we do not attempt to use this approach 
in combination with the measured EW of the lines to estimate the 
metallicity of the gas and constrain the hydrogen column density. 

In our analysis we only considered absorption systems with unsat- 
urated OVI, OVII, OVIII and NelX lines with EW > 20 km 
i.e. strong enough to be detected. For the OVI we have also con- 
sidered weaker lines with EW > 5 km potentially observable 
in the UV spectra. The density and temperature of the gas can be 
found by minimizing the difference between EWp {X^ ) /EWp {Y-' ) 
and EWm{X'') / EWm{Y-' ), and considering at least two indepen- 
dent EW ratios. Figure[TT]illustrates the outcome of this minimiza- 
tion procedure. Each curve represents the locus of the minima for 

the function [|^{|^ - eZI(y?) (fe' ^g)] in the phase space, 

given the measured ratio EWm{X^)/ EWm{Y-'). The two curves 
refer to the OVIII/OVII and NelX/OVII ion ratios. The estimated 
value for the gas density and temperature is found at the intercep- 
tion between the two curves. The cross indicated the "true" density 
and temperature of the absorbing gas. As already discussed in Sec- 
tion 3.1 these values of pg and Tg are somewhat ill-defined since the 
physical state of the gas can vary across the absorbing region and the 
effective values of its density and temperature depend on ion con- 
sidered and on the estimator adopted to measure its opacity. This 
ambiguity can be regarded as an intrinsic uncertainty in the density 
and temperature of the absorbing gas. In Section 3.1 we have used 
two different estimators (Eq. |2] and [3]l to quantify this uncertainty 
and found that the typical difference in both temperature and density 
found when using either estimator is 15-20 % represented by the size 
of the cross in Fig.E] To further assess the amplitude of this scatter 
we have used yet another estimator in which the temperature and the 
density of the gas are measured in correspondence of the minimum 
of the absorption line. Using the full battery of the three estimators 
confirms that intrinsic uncertainty in the density and the temperature 
of the absorbing gas is 15-20 %. 

The results presented in this Section refer to the mock absorbers 
for which we can measure the OVII/OVIII and NelX/OVII EW ra- 
tios. Solutions to the minimization procedure are searched for in 
the phase-space region with lO"*-^ < Tg < lO'' K and lO"'^ < 
Pg < 10^^ cm" , corresponding to over-densities in the range 
5 < (5g < 5000, i.e. slightly larger than the canonical WHIM range. 
In ~ 70 % of the cases we have been able to estimate the physical 
state of the absorbers. For the remaining ~ 30 % of the absorbers, 
we fail to find a common solution, i.e. the two curves do not cross 
in the explored density-temperature range. We have checked that if 
we include the third, non-independent curve relative to EW ratio of 
NelX/OVIII, the resulting additional crossing points are close to the 
original one. 

To assess the accuracy of our procedures in the 70 % successful 
cases, we have compared the pg and Tg estimated through the min- 
imization procedure with the true values. Random and systematic 
errors can be appreciated from scatterplots like those in Figs.ll4land 
[15] that refer to the mock absorption systems featuring detectable 
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Fig. 11. — The locii of the minima for the function 

[ flWM(yJ) ~ fiWp(yj) ] "^'^ (Ps^Tg) plane. The two curves 
refer to a mock absorber at z = 0.418 for which we have detected the OVll, 
OVlll and NelX hnes. The cross indicates the expected values. 

OVII, OVIII and NelX lines. The estimated gas density turned out 
to be systematically smaller than the true one by about 35 %, sim- 
ilar in size to the random error of 50 % that we identify with 
the scatter around the best fitting line in the plots. The temperature 
is overestimated by about 35 %, with a random error of the same 
amplitude. Using OVI instead of NelX leads to similar, although 
somewhat smaller, systematic and random errors. 

These errors may originate either from the presence of strong in- 
homogeneities within the absorbing gas or from the breakdown of 
the optically thin approximation. To check the validity of the first 
hypothesis we have examined the behavior of the temperature and 
the density within the redshift range of the lines corresponding to the 
WHIM absorbers. Fig.[T2]is analogous to Fig. |2]and shows a typical 
unsaturated mock absorber. The temperature and the density of the 
gas are remarkably constant across the redshift range of the OVII and 
OVIII lines. This result is at variance with that found by Kawahara 
et al. (2006). The authors found large deviation from homogeneities 
in the absorbing gas, as clearly illustrated in Fig. 15 of their paper. 
This discrepancy can be understood by noticing that Kawahara et al. 
(2006) absorbers are typically found in regions of higher density and 
temperature than ours, i.e. within environments in which density and 
temperature vary change significantly across a ~ Mpc scale. Such 
a difference probably derives from the fact that the OVII ionization 
fraction and the overall gas metallicity in Kawahara et al. (2006) 
WHIM model are smaller than in our case. 

To test whether the systematic errors in the estimate of pg and Tg 

can be attributed to departures from the optically-thin approximation 

we have computed fractional difference between the measured and 

the predicted EW ratio computed using the true values of density 

and femneramre A - \ EWp{N^lX)EWM{Ovn) _ ,] ^, , , 
ana temperature, z:i - yEWp(ovii)EWM(Neix) • ine DieaK- 

down of the optically thin approximation leads to overestimate the 
equivalent width of the line, EWp > EWm- This mismatch is 
more serious for OVII, which is the strongest line, and should in- 
crease with the strength of the line. As a consequence we expect 
a negative A with an amplitude that increases with EWm{OVII). 
This effect can indeed be appreciated from Fig.[T3]in which, super- 
imposed to the scatter, there is a clear trend for A to decrease with 
EWm{OVII). The scatter is particularly large for weak lines. It is 
contributed by the fact that the validity of the approximation does not 
solely depends on the EW of the OVII line and, to a lesser extent, 
by the presence of inhomogeneities in the absorbing gas. Although 
we show only the case of NelX and OVII, we have checked that the 
same trend is found for the other ratios that use the OVIII and OVI 



Studying the WHIM with GRBs 



13 




1 

0.8 
0.6 
0.4 

7 
6.5 

6 
5.5 

5 
4.5 

-3 
-4 
-5 
-6 
-7 





- , , , , 1 


1 




1 


2 


3 


,1111, 
'- 


1 






1 


2 


3 


.1111,1111,1111. 


: , , , , 1 


1 





1 2 
l.o.s. distonce [Mpc h"^] 



Fig. 12. — Example of a typical unsaturated WHIM absorber in our mock 
spectra from model Bl. From Top to Bottom: Transmitted flux for OVll, 
dVlll, temperature and density of the gas along the LOS distance. 

lines. We note that the ~ 30 % of the absorbers for which the mini- 
mization fails to find a (pg, Tg) solution typically have a large A. 

The existence of an anti-correlation between A and 
EWm{OVII) suggests a practical, although approximate 
way to correct for systematic errors since one can use the A - 
EWm{OV 1 1) relation in reverse to correct the predicted EW ratio 
knowing the EW of the OVII line. 

To first approximation, this correction should not depend on the 
WHIM model since for a given gas density, temperature and relative 
metal abundance, the predicted EW ratio depends only on the ion- 
ization model. Therefore, under the hypothesis of hybrid, photoion- 
ization plus collisional ionization, and assuming the same UV and 
X background we expect that the mean A - EWm{OV 1 1) relation 
should not depend on the WHIM model. To check this hypothesis 
we have repeated the same analysis using model B2 instead of Bl 
and found that, at least in the range of EWs we are interested in, 
the linear fit to the mean A — EWm (OVII) relation is consistent, 
within the errors, with the one shown in Fig. 1131 

Figs. [14] and 1 151 compare the logarithm of the density and tem- 
peratures predicted from the OVII, OVIII and NelX EW line ratios 
to their true values after applying the statistical correction described 
above. The correction is very effective in reducing systematic but 
also random errors. Systematic errors on the density and the temper- 
ature decrease to ^ 17 % and ~ 10 %, respectively and are compa- 
rable, to random errors of ~ 23 % and ~ 25 %. Similar results were 
obtained for absorbers featuring the OVI rather than the NelX line. 
In this case, however, systematic errors are further reduced to a few 
percent level. 

To summarize, our analysis shows that one can measure the effec- 
tive temperature and density of an absorption system by measuring 
the EW ratio of at least three lines. However, the estimate is affected 
by systematic errors that are mainly contributed by the breakdown of 
the optically-thin approximation. Nonetheless, once the EW of the 
OVII line is measured accurately, unbiased estimates can be obtained 
through simple and robust corrections that depend on the ionization 
model adopted. Applying this technique to absorption systems fea- 
turing the NelX line detected in the GRB afterglow spectra measured 
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lines with EW > 20 km are shown with black, filled dots. 

by EDGE/XENIA, one expects to be able to measure the density and 
temperature of about 16 absorbers per year with ~ 25 % random 
errors. This estimate is obtained by multiplying the number of de- 
tectable absorption systems listed in Table 4 by the the 70 % success 
rate of the estimation technique. 

We note that this estimate requires a-priori knowledge of the Ne/0 
relative abundance. This needs not to be true when we focus on those 
absorbers featuring the OVI, OVII and OVIII lines. Unfortunately, 
because of the low OVI column density, the OVI ls-2p line is very 
weak and one expects to detect only ~ 2 such absorbers per year, 
as shown in Table 4. Once more, to increase the number of OVI 
detections one should resort to simultaneous UV spectroscopy. 

Finally, we point out that we have only considered Neon and Oxy- 
gen ions. Using all other ions listed in Table 4 would significantly in- 
crease the the number of absorbers with at least three lines for which 
one can estimate the temperature and the density of the gas. 

5.2. Measuring the WHIM cosmological density 

Another interesting question is how to estimate the cosmological 
mass density of the WHIM from the line EWs measured in the GRB 
afterglows spectra. As pointed out by Richter et al. (2008) this 
estimate can be obtained in two steps. First, following the procedure 
outlined in the previous Section, one uses the measured EWs of the 
ion lines and their ratios to determine the ion fractions and, from 
these, the total gas column density for a given metallicity. Then, 
one integrates the gas column density along LOS to different X-ray 
sources in order to minimize statistical errors and cosmic variance. 
Following this procedure, the cosmological density of the WHIM in 
units of the critical density pc SIwhim, can be estimated from all 
absorption systems featuring the OVII line as: 



V f Ah.'\ ( N(OVII) \ 
^''^ Uo A,, \X{OVII))^, 



(7) 



where p — 1.3 is the mean molecular weight and mn is the mass 
of the hydrogen atom. The index i denotes the absorption systems 
along the LOS j. Each absorber is characterized by the OVII col- 
umn density Novii, that can be evaluated from the line EW, the 
ion fraction, Xovii, which can be evaluated from the EW ratios, 
and the metallicity {Ao /Ah) that needs to be assumed a priori. The 
quantity dc{zj) denotes the comoving distance corresponding to the 
redshift range in which absorption systems can be detected along the 
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Fig. 14. — Scatterplot of the predicted gas density, pp vs. measured gas 
density pM after con'ection for systematic eiTors. Filled dots refer to all ab- 
sorbers with EW (OV 1 1) > 20 km s^^ and filled, green squares to ab- 
sorbers which are saturated in any of the lines. The dashed red line represents 
the best fit to the unsaturated points. 
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Fig. 15. — Scatterplot of the predicted gas temperature, Tp vs. measured 
gas temperature Tm after correction for systematic errors. Filled dots refer 
to all absorbers with EW{OVII) > 20 km and filled, green squares 
to absorbers which are saturated in any of the lines. The dashed red line 
represents the best fit to the unsaturated points. 

line of sight j and depends on the assumed cosmology. It is straight- 
forward to generalize Eq. |7]to account for all remaining absorption 
systems with featuring lines different from OVll 1 which, however, 
are very few. 

Nicastro et al. (2005) used an estimator similar to that of Eq. [7] 
to measure SIwhim from the putative absorption systems along the 



LOS to Mrk 421. They assumed that the {Ao /Ah) ratio is constant 
and that errors are dominated by Poisson noise from low number 
statistics. Under the same hypotheses, one would expect to estimate 
f^WHiM with ~ 10 % accuracy using the Naba ~ 50 absorption sys- 
tems that EDGE/XENIA is expected to detect in the GRB afterglows 
spectra during a 3-year observational campaign. 

This error estimate is clearly too optimistic since it does not ac- 
count for the uncertainties in the estimated ion fraction and ignores 
both the cosmic variance and the fact that different absorbers may 
have different metallicity. In an attempt to obtain a more realistic er- 
ror estimate we have measured J^whim using all absorption systems 
listed in Table 4 through the following, simple estimator: 



WHIM = 



1 



Pc 



(8) 



p® represents the effective gas density of the absorber estimated from 
the EW ratios using the technique described in the previous Sec- 
tion and for which we have also estimated the associated errors. 
(I^whim) represents the average comoving depth of the WHIM ab- 
sorbers. This quantity has been computed using all mock absorption 
systems with a detectable OVlI line and represents the comoving dis- 
tance corresponding to the redshift interval [z, 2] in which the trans- 
mitted OVlI flux is below 95 %. Clearly, (Dwhim) is expected to 
depend on the WHIM model adopted. For model 51 we found that 
(fwHitvi) = 1.4 Mpc with a scatter around the mean of 0.6 Mpc 
that quantifies the random error Uncertainties on SIwhim have been 
obtained by adding in quadrature the errors in the measured and 
the scatter in Dwhim and and by assuming that the total uncertainty 
decreases as \/N abs- We neglect the cosmic variance. Under these 
hypotheses the errors on JIvvhim is ~ 15 %, not too different from 
the simple Poisson estimate. 

Applying Eq. [8] to the mock absorption systems we obtain 
f^wHiM = 0.014 ± 0.002, slightly smaller that the expected value 
f^wHiM = 0.017. The 1.5 ct difference reflects the fact that the 
distribution of Dwhim is not Gaussian but is skewed toward large 
values. 

This exercise suggests that the analysis of X-ray GRB spectra 
taken in 3 years by a satellite like EDGE/XENIA should allow to 
measure Owhim with 15 % accuracy. However, we stress that this 
estimate strongly depends on the WHIM models adopted, which is 
used either to determine the metallicity of the absorbers, [Ao /Ah)- 
in Eq.|7]or their average depth, (Dwhim, ) in Eq.[8] 

6. DISCUSSION AND CONCLUSIONS 

In this work we have considered different numerical models of the 
WHIM consistent with currently available observational constraints 
to assess the possibility of detecting the signature of the missing 
baryons in the X-ray spectra of distant sources and to study their 
physical properties. In particular we focused on the absorption spec- 
tra of distant GRB afterglows that could be taken by the Wide Field 
Spectrometer on board of the recently proposed EDGE and XENIA 
missions. Both satellite are capable of fast repointing, which make 
them suitable for observing the GRB afterglows during their early, 
bright phase. We note that although we have only considered ab- 
sorption spectroscopy the WFS is especially suitable for studying 
the WHIM in emission, as we will show in a future work. The main 
results of this study can be summarized as follows: 

(i) All WHIM models explored provide similar predictions for 
the cumulative distribution of OVI line EW but not for OVII and 
OVIII. This fact reflects the differences in the metallicity prescrip- 
tions among models, which become more evident in regions of en- 
hanced density and temperature in which the OVII and OVIII ions 
are more common than OVI. As a consequence the OVI lines are not 
efficient signposts for the bulk of the WHIM which, instead, is best 
traced by OVII and OVIII, as already shown by Chen et al. (2003). 

(ii) GRBs are better X-ray background sources than AGN in many 
aspects. First of all, for a reference fluence in the [0.3-2] keV band 
of lO"*^ erg cm"^ requked to detect the WHIM, the typical GRB 
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distance is larger than that of AGN, hence increasing the probability 
of WHIM detection. Second of all, GRBs constitute a replenishable 
reservoir of sources, contrary to AGN. As a consequence, the total 
number of GRBs above the reference fluence that EDGE/XENIA is 
expected to detect in a year is already a factor of two larger that the 
total number of available AGN with z > 0.3 with the further advan- 
tage that the photons can be collected in half observational time. 

(Hi) A satellite like EDGE/XENIA is expected to detect 25-60 
WHIM absorbers per year along the LOS to ~ 30 GRBs through 
the OVII line. The statistical significance of the detections is 4.5 a 
which guarantees that the probability of spurious line contamination 
is less than 1 %. The lower bound corresponds to the most conserva- 
tive WHIM model considered here (dubbed A/) which, however, we 
regard as less realistic since it assumes a deterministic metallicity- 
density relation. The upper value derives from a more physically 
sound WHIM model (-82) that accounts for the scatter in the rela- 
tion. A more realistic upper bound should also account for the pos- 
sibility of detecting WHIM absorbers beyond z — 1, which we have 
ignored in our analysis. Finally, with a goal energy resolution of leV 
instead of 3 eV (baseline) the number of expected detections would 
further increase by about 40 %. 

(iv) A significant fraction of the WHIM absorbers detected 
through the OVII line comes with an additional, detectable OVIII 
line. The total number of expected WHIM absorption systems fea- 
turing both the OVII and OVIII line is about 35 per year when we 
consider model B2. The simultaneous detection of both lines in- 
creases the significance of the detection to 5 cr. A large fraction of 
these systems have a third, detectable ion line, typically of a different 
element, that further increases the significance of the detection. 

(v) Measuring the EW of at least three lines allows to characterize 
the thermal properties of the absorber. Under the condition of hybrid 
collisional and photoionization equilibrium adopted in our WHIM 
models and assuming an optically thin medium, we were able to de- 
termine the effective density and temperature of the gas in ^ 70 % of 
the line absorption systems featuring a NelX or OVI line in addition 
to OVII and OVIII. These estimates, however, are affected by sys- 
tematic errors that derive from the breakdown of the optically thin 
approximation. A simple, model independent correction can be ap- 
plied that allows to remove this systematic effect. As a result one 
should be able to estimate the logarithm of the temperature and den- 
sity of about 20 WHIM absorption systems per year with ^ 25 % 
random errors. 

(vi) Integrating the gas density along the lines of sight to the mul- 
tiple line absorbers that a mission like EDGE/XENIA is expected to 
detect during the first three year of the mission, would allow to esti- 
mate the cosmological mass density of the WHIM with an accuracy 
of ~ 15 %. This estimate, however, heavily relies on the assumed 
gas metallicity which currently constitutes the main uncertainties of 
the WHIM models. 

As soon as they will become available. X-ray observations will 
be used to validate the WHIM models and therefore to constrain the 
poorly known metal enrichment mechanisms. Nevertheless, with no 
rigorous metallicity model available, the abundance of the heavy ions 
in the WHIM can only be measured by observing the Ly-a hydrogen 
line associated to the X-ray absorbers. The Cosmic Origin Spectro- 
graph will allow to search for broad Ly-a absorbers which are likely 
associated to the WHIM along the LOS to bright AGN. For tran- 
sient sources, like GRBs, a similar study would require the possibil- 
ity of performing simultaneous X-ray and UV spectroscopy which is 
currently not envisaged in the EDGE or XENIA mission concepts. 
Yet, the success of a combined UV and X-ray spectroscopy strategy 
strongly depends on how well broad Ly-a absorbers trace the WHIM 
gas. Current datasets provide no clue on their mutual relation. On 
the theoretical side Chen et al. (2003) found weak correlation be- 
tween HI and both OVII and OVIII that, as we have stressed already, 
trace the bulk of the WHIM. However, a dedicated theoretical study 
must be performed to assess how well a combined UV and X-ray 
spectroscopy will measure the metal content of the WHIM gas. 

Although we did not mention it explicitly, the Wide Field Spec- 
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Fig. 16. — Bottom panel. Zoo med plot in a narrower energy range of the 
mock absorption spectrum in Fig. 1101 centered on the OVII - OVlll absorp- 
tion systems at 2 = 0.069. Top panel: Emission spectrum in a 4 armin^ area 
along the same direction, taken in a 1 Ms observation with EDGE/XENIA 
Wide Field Spectrograph. The OVII triplet and the OVIII resonant line asso- 
ciated to the absorption lines are clearly seen along with the Galactic OVII 
triplet. 



trometer on board of EDGE or XENIA offers the possibility of com- 
bining emission and absorption analyses. This is an important asset 
of the EDGE/XENIA mission concepts since detecting the OVII res- 
onant line in absorption and the triplet in emission allows to constrain 
the metal abundance of the WHIM. As in the case of the combined 
UV and X-ray spectroscopy, the success of this alternative strategy 
depends on how well WHIM absorbers can trace the emitting gas 
and on our ability to remove or model the contribution of the differ- 
ent backgrounds/foregrounds. The combined study of the WHIM in 
absorption and emission is the subject of future work. Here we antic- 
ipate that on average 30 % of the line systems detected in absorption 
can also be seen in emission in deep observations. Figure [16] illus- 
trates one of these cases. The lower panel shows a zoom in a narrow 
energy interval of the mock spectrum displayed in Fig.[TOl centered 
on the OVII and OVIII lines of the absorption system at z — 0.069. 
The upper panel shows the emission spectrum in a 4 armin^ area 
taken along the same LOS when the GRB is off in a 1 Ms long ob- 
servation. The OVIII and OVII lines are detected both in absorption 
and in emission. All components of the OVII triplet are clearly seen 
and can be well separated from the Galactic component. The com- 
bined emission/absorption spectroscopy not only constitutes a possi- 
ble way to constrain the WHIM metallicity and hence to estimate its 
cosmological density, but also suggests an observational strategy for 
the WHIM which consists in pre-selecting target fields in which the 
WHIM absorption lines have been detected in the X-ray spectrum of 
a GRB for deeper, follow-up observations once the GRB has faded 
away. 

The possibility of studying the WHIM in emission is crucial to 
investigate its spatial distribution, since the pencil-beam-like sam- 
pling provided by absorption spectroscopy is too sparse to trace its 
large scale, filamentary structure that we observe in the simulations. 
Yet, cross-correlation of WHIM absorbers and next generation all- 
sky galaxy surveys out to z = 2 like ADEPT or EUCLID may shed 
light on the WHIM-galaxy connection, very much like in the recent 
correlation analyses of weak Ly-a absorbers and Hl-rich galaxies 
performed by Ryan- Weber (2006) and Pierleoiri et al. (2008) and 
that of OVI absorbers and galaxies (Ganguly et al. 2008). These 
studies have helped in clarifying the relation between the diffuse 
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baryons in the intergalacdc medium and tlie atoms locked inside viri- 
alized structures. 

Finally, we would like to stress that although X-ray spectroscopy 
constitutes the best way to study the missing baryons, alternative 
strategies like the cross correlation of the kinetic Sunayev-Zel'dovich 
signal near clusters and groups and the all-sky E-mode polarization 
in the context of upcoming CMB experiments like Planck, ACT, SPT 
or APEX will provide useful, independent constraints on the pres- 
ence and abundance of missing baryons in the Universe (Hemandez- 
Monteagudo & Sunyaev 2008). 
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